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Comment on "A simple one-dimensional model of 
heat conduction which obeys Fourier's law" 

In a recent letter Garrido et al Q consider heat con- 
duction in a model of hard point particles of alternating 
masses on a one-dimensional line. Based on a number 
of numerical results, the authors claim that this momen- 
tum conserving model exhibits Fourier's law. We first 
comment on the apparent contradiction with an earlier 
result of Prosen and Campbell || (PC). We then point 
out certain inconsistencies in their results and disagree- 
ments with our own results. 

The authors have first measured the system size de- 
pendence of the mean current (J) = (N _1 ^ miuf/2) 
where mi, xi and u\ denote the mass, position, veloc- 
ity of the /th particle and N is the number of parti- 
cles . As they correctly point out, it is not possible 
to make definite conclusions from this simulation data, 
since the asymptotic regime may not have been attained. 
Next the authors compute the current-current correlation 
C(t) = N(J(t)J(0)) and find a decay C*(t) ~ i -1 - 3 , which 
is sufficiently fast to give a finite Kubo conductivity k. 
This would seem to contradict the exact result of PC on 
infinite thermal conductivity in momentum conserving 
systems. Their proof applies to this model. This ap- 
parent contradiction can be explained by the fact that, 
in their simulation, Garrido et al work in an ensemble 
with total momentum set to zero and in this case the PC 
proof does not predict anything. As has been pointed 
out by Bonetto et al 0], the correct expression for the 
Kubo formula requires one to use the connected part of 
the correlation function if one is working in the canonical 
ensemble [0. Alternatively one may fix the momentum 
to be zero and work with the usual correlation function 
as Garrido et al have done. Thus the work of PC does 
not really prove divergence of k. 

However some aspects of the paper appear to be unsat- 
isfactory and need some explanation. Firstly the linear 
temperature profiles obtained in the paper are inconsis- 
tent with the finding of finite k. The temperature de- 
pendence of k can be scaled out from the Kubo formula, 
yielding a T 1 / 2 dependence (given k is finite). This fol- 
lows from the fact that the correlation Cr(t) at temper- 
ature T has the scaling form C T (i) = T*Ci(T l IH). The 
T 1 / 2 dependence of k also follows from simple kinetic the- 
ory arguments. A temperature dependent k at once leads 
to nonlinear temperature profiles. Infact in our study of 
the same model || we clearly see the expected nonlin- 
ear profiles for similar system sizes. One possible reason 
for the difference could be that Garrido et al use deter- 
mistic heat baths unlike the stochastic heat baths used 
by us. It is not clear how well such deterministic baths 
simulate true thermal sources. Another possibile source 
of error is the way Garrido et al define local tempera- 
ture, namely, by measuring the mean kinetic energy and 
mean position of each particle. In one dimensions posi- 
tion fluctuations are large (~ y/N) and a more correct 



procedure is the one used by us: to define the local num- 
ber and energy densities as n(x, t) — 5{x ~ x l)) an d 
e(x, t) — (J2i( m i u2 /2)H X ~ x l)) respectively, and then 
define the local temperature as T(x) = 2e(x)/n(x). 

Secondly, from our own simulations, we are unable to 
verify the results of fO]. The authors have computed C(t) 
and also the onsite correlator c(t) = (ji(t)ji(0)) . They 
find that for the unequal mass case, C(t) and c(t) have 
roughly the same long-time decay ~ 1/t 1 ' 3 while for the 
equal mass case c(t) ~ 1/t 3 . Our results are summarized 
in Fig. (0) and the important differences with |jj are: 

(i) We do not find any evidence for the decay C(t) ~ 
1/t 1 - 3 . Infact the behaviour we find is consistent with 
the decay J - 1/N°- S3 found in [§. However we do not 
wish to make a strong case for the number 0.83 though 
it is rather striking that we get the same number from 
two very different approaches. 

(ii) The behaviour of c(t) seems to be very different 
from that of C(t) contrary to what is claimed in The 
authors comment that c(i) has better averaging proper- 
ties and, because it shows approximately the same decay, 
this confirms the behaviour seen for C(i). But is there 
any reason to expect that for the unequal mass case C(t) 
and c(t) will have similar decay laws ? Infact for the 
equal mass case we know that C(t) is a constant (since J 
is a constant of motion) and therefore the behaviour of 
C(i) and c(i) are drastically different. 

(iii) The equal mass case is nonergodic since there are a 
macroscopic number of conservation laws. Thus time av- 
erages from simulations are initial-condition-dependent. 
We can verify this in our simulations and also find that 
making the masses slightly unequal restores ergodicity. 
Thus it is hard to understand how the decay c(i) ~ 1/t 3 
is obtained by |Q . We note that the paper by Jepsen B 
(quoted by Garrido at al) only gives (vi(t)vi(O)) ~ 1/t . 
It is not clear how this leads to the same prediction 
for the asymptotic behaviour of c(t) which is more like 
(vf (t)v 3 (0)) . Besides, the calulation of Jepsen is not done 
in the zero-momentum ensemble. 

In our simulations time averages were performed over 
10 9 — 10 10 collisions. As checks on our simulations we 
found that C(0) and c(0) agree with their exactly known 
values and also that C (t) and c(t) satisfy the exact scaling 
relations mentioned above. In the equal mass case we got 
a constant C(t), as expected. 

In conclusion we have shown that there does not seem 
to be any hard evidence to prove validity of Fourier's 
law in this system. We think it is as hard to conclude 
anything (about finiteness of k ) from correlation function 
data as it is from the Jversus N data. 

I thank Onuttom Narayan for illuminating discussions 
and for critically reading the manuscript. 
Abhishck Dhar 
Raman Research Institute 
Bangalore 560080 
India. 
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FIG. 1. Plot ol C(t) for system sizes N = 100 
(dot-dashed), 200 , 400 and 800(solid) [ T = 1, mi = 1, 
m 2 = 2]. Inset shows c(t) for N = 100 (dots), 200 and 400 
(dot-dashed) The straight line corresponds to the power-law 
decay ~ l/t ' 83 . 
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